Quantum Transport of Massless Dirac Fermions 
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Motivated by recent graphene transport experiments, we have undertaken a numerical study of 
the conductivity of disordered two-dimensional massless Dirac fermions. Our results reveal distinct 
differences between the cases of short-range and Coulomb randomly distributed scatterers. We 
speculate that this behavior is related to the Boltzmann transport theory prediction of dirty-limit 
behavior for Coulomb scatterers. 
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Introduction — Graphene can be described at low- 
energies by a four component massless Dirac-fermion 
(MDF) modelfT] that has long attracted theoretical at- 
tention because of appealing properties including chi- 
ral anomalies0J^, 0, randomness induced quantum 
criticalit yla |7|, unusual electron-electron interaction 
effects [ol-lioi. relevance to high Tc superconductors 
and various unusual transport propertieslEBBr 
HI, H E, [11 E H [13, [ii, [3, [3 • The recent experi- 
mental realization of single-layer graphene sheets [2l| has 
made it possible to confirm a number of theoretical pre- 
dictions, including unusual quantum Hall effects 22, |23|, 
and has also revealed some surprises. The main find- 
ings can be summarized as follows: i) graphene's conduc- 
tivity a never falls below a minimum value (cr""") cor- 
responding (approximately) to a conductance quantum 
(e^//i) per channel, in spite of predictions [gI. [tI. [ifl. [l^ [l6l| 
based on the self-consistent Born approximation (SCBA) 
that fj™™ = (I/tt) e'^/h for a MDF channel, and predic- 
tions that localization occurs [l^ [13, [l3l when interval- 
ley scattering is significant; ii) in gate-doped graphene 
a increases linearly with carrier density n away from 
the charge neutrality (Dirac) point, implying a constant 
mobility /i = i7/ne[22. l23j and not the constant con- 
ductivity usually predicted for the Boltzmann transport 
14| . Although these surprises have inspired a 



regime,^ 

number of theoretical studies [3, 17, [l^, the source of 
the discrepancies between experiment and theory has not 
yet been conclusively identified. 

We have recently pointed out that the linear depen- 
dence of conductivity on carrier density in graphene can 
be explained in the framework of Boltzmann transport 
theory by assuming Coulomb scatterers [2^ rather than 
the short-range scatterers assumed, mainly as a practical 
simplification, in most theoretical work. The golden-rule 
scattering rate for Coulomb scatterers in a MDF model 
diverges in the zero energy (Dirac point) limit, whereas 
it vanishes for short-range scatterers. This property sug- 
gests the possibility of a qualitative difference between 
these two disorder models. Since Boltzmann theory is 
not applicable in the vicinity of the Dirac point, how- 
ever, a fully quantum mechanical approach is required to 
address the minimal conductivity. In this Letter we re- 



port on a finite-size Kubo formula analysis used to study 
the quantum transport of MDFs. Over the range of sys- 
tem sizes that we can describe, we find that ct™'" is a few 
times larger for Coulomb scatterers than for short-range 
scatterers and ^ e'^/h. 

A MDF model describes graphene transport only when 
intervalley scattering is unimportant. We argue that the 
MDF model likely does apply at accessible experimental 
temperatures to graphene systems near the Dirac point, 
since the intervalley scattering length obtained within 
the Born approximation diverges in that limit. The tem- 
perature below which intervalley scattering is important 
should increase away from the Dirac point ^25sj. 
Massless Dirac Fermion Model — Graphene's honeycomb 
lattice has two atoms per unit cell on sites labeled A and 
B. The low-energy band structure consists of Dirac cones 
located at the two inequivalent Brillouin zone corners K 
and K': 
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and iJx' — fiver*' ■ k, where v is the graphene Fermi ve- 
locity and the Pauli matrices cr act on the sublattice 
degrees of freedom. For each wavevector k, Eq.lJT]) has 
two eigenstates |k,±) = (|k,A) ± e"^|k, B))/\/2, where 
4> = iavT^ {ky / kx) , and eigenenergies E-^^ — ±hv\\i\. 
When intervalley scattering is neglected the K and K' 
valleys contribute independently to the conductivity. 
Boltzmann theory for doped graphene — We start by 
briefiy reviewing Boltzmann transport theory applied 
to graphene since this consideration motivates Colomb 
scatterer models. The Boltzmann conductivity CTq — 
(e^//i) {2EpTQ/h) = {e? jh) 2kp£ is proportional to the 
transport relaxation time Tq. For short-range scatter- 
ers the Born approximation gives H/tq = 2itV^pf — 
(n^u^/hv^) \Ef\ [ii, where V{r) = u^f (^(r-R/) 



is the disorder potential, is the density of scatterers, 
and pp the density of states at the Fermi level. When 
the range of the impurity potential is much longer than 
the lattice s pac ing of graphene, intervalley scattering is 
weak[illil,|i7|. 

Note that ctq is independent of the carrier density n. 
Experiment, on the other hand, finds that the mobil- 
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tion, and |n) denotes an eigenstate of the Dirac equation, 
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FIG. 1: (Color online) Dirac fermion conductivities for (a) 
short-range scatterers and (b) screened Coulomb scatterers. 
The inset of (a) compares the densities of states for short- 
range and Coulomb cases, (c) shows Kubo formula conductiv- 
ities at the Dirac point as a function of ri for Coulomb scatter- 
ers for three different system sizes. The inset of (c) shows the 
size dependence of the conductivity for unscreened Coulomb 
(open circles) and screened Coulomb (boxes) cases, (d) Distri- 
bution function of the conductivity and resistivity (inset) at 
the Dirac point for Coulomb scatterers. In the main panel, 
the solid (blue) line is for L — lALmin and the dotted (red) 
line is for L — 2Lmm, where Lmin is the minimum system size 
considered. 



ity n = a/ne in graphene is nearly constant except at 
very low-densities. One plausible explanation for this 
behavior is that Dirac fermion scattering is dominated 
by Coulomb scattering from ionized impurities near the 
graphene plane, V{t) = ^2?' e^/^lr — R./|. Using Fermi's 
golden rule and approximating the screened Coulomb 
interaction by'ii'l J7sc(q) = (27re^)/e(g -I- -iaghp) ~ 
(HvTr)/ (2kp), the Boltzmann conductivity for Coulomb 
interactions is erg ~ {Ae'^ /h) {n/rii) 32/7r, proportional to 
density in agreement with experiment. Here Ug = e'^/ehv 
is the effective fine structure constant used to character- 
ize the ratio of Coulomb interaction and band energy 
scales in graphene. {ag ~ 3 in vacuum and ~ 1 when the 
graphene sheet is placed on a Si02 dielectric substrate.) 
Note that the Boltzmann conductivity for Coulomb scat- 
terers vanishes as n — > 0, contradicting experiment [2^ . 

MDF model finite-size Kubo Formula - Our numerical 
results, obtained by evaluating the finite-size Kubo for- 
mula 



ihe' ^ f{E, 

n.n' 



/(£'„/) (■n\vx\n'){n'\vx\n) 

En' En - En' + IT] 



(2) 



ihva ■ V + V{v) 



^ Eip 



(3) 



are summarized in Fig.l. Here v = vcr is the Dirac 
fermion velocity operator, f{E) is the Fermi-Dirac func- 



which we solve using a large momentum-space cutoff A. 
The disorder potential momentum-space matrix elements 
are 

1 ^" 

(ka|y|kV) = -^Y.^(^~ e'^^-^'^-"^', (4) 

1=1 

where U{q) is constant for delta- impurities and given by 
C/sc(q) for the screened Coulomb scattering case. The 
scattering center locations R/ and potential signs were 
chosen at random. The results shown in Fig.l are typ- 
ically averaged over 10"* disorder realizations. We es- 
timate the bulk conductivity by evaluating Eq.([2]) at a 
large number of r;- values. The Kubo conductivity a van- 
ishes for both small and large rj but there is an inter- 
mediate region where the dependence of cr on 77 is rel- 
atively weak. We use the maximum of a vs. rj to es- 
timate the conductivity at a given system size L. For 
metals, including doped graphene, physical arguments 
suggest that 77 ^ ^/Tl where Tl is the escape time 
from the system studied numerically, ft follows [26j that 
77 ^ (A-E) ~ gT 5E where the Thouless energy {IS.E) 
is the geometric mean of the eigenvalue difference be- 
tween periodic and antiperiodic boundary conditions, 
6E — 1/L'^pp is the level spacing at the Fermi level, 
and grp = {AE)/SE is the Thouless number. 

The reliability of the finite-size Kubo formula method 
described above is solidly established for diffusive metal- 
lic systems [ill. For conventional non-relativistic two- 
dimensional electron systems with parabolic dispersion, 
numerical conductivities calculated in this way agree 
closely with the theoretical expectation of a Drude 
conductivity (Tq — {e^ lh)kF(. with a finite-size weak- 
localization correction ^23|. Subtleties, related to the rj 
dependence of the numerical estimate, do however arise 
when this method is applied to insulators. The Dirac 
point of the MDF model is an intermediate case and re- 
lated uncertainties apply to our numerical estimate of 
ct'°'". The property that the Dirac point density of states 
is finite in the presence of disorder, as illustrated in the 
inset of Fig. 1(a), may help vahdate the Kubo approach. 
We find that the dependence of cr on at the Dirac point, 
illustrated in Fig. 1(c) is similar to that in a metal[2^. 
indicating delocalized states, cr is a maximum when rj 
corresponds to ~ 5E at each system size as expected for 
a ^ e^ /h. We have also evaluated the Thouless conduc- 
tivity as a consistency check. A large number of numer- 
ical studies have demonstrated that the Thouless con- 
ductivity estimate, although perhaps not as accurate in 
the diffusive metal limit (cr >> e^//i), is more univer- 
sally applicable. It may be used for both delocalized and 
strongly localized states[2^, and should therefore be re- 
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liable at the Dirac point. We find consistency between 
the two a estimates. 

Fig.l compares Kubo conductivity estimates for (a) 
the short-range scatterer and (b) the screened Coulomb 
scatterer cases. For the short-range disorder potential 
case, the density dependence of the conductivity is non- 
linear, approaching the constant Boltzmann conductivity 
for \Ef \ ^ h/r. The estimated value of a{E = 0) is close 
to the SCBA value (l /Tr)e'^/ h predicted in earher theo- 
retical studies 14, For the screened Coulomb 
scatterer case, on the other hand, the conductivity a in- 
creases linearly with increasing density |n| as Boltzmann 
theory predicts. At the Dirac point, however, the con- 
ductivity remains finite with the minimum value ~ e^//i, 
which is few times larger than the SCBA value for the 
short-range model. To illustrate the dependence on car- 
rier density, we have smoothed the curves in Fig. 1(a) and 
(b) by averaging over Fermi energy interval containing 
typically 1-30 levels, and over boundary conditions, in 
addition to over approximately 10^ disorder realizations. 

If we neglect intervalley scattering and account for 
graphene's spin and valley degeneracies, MDF proper- 
ties can be compared with graphene experimental results. 
Both the linear dependence of the conductivity on density 
and the shift of a™™ suggest that long-range scattering 
similar to that produced by ionized impurities is present 
in experimental samples. We note that ionized impuri- 
ties in the substrate or at samples edges that are sep- 
arated from the conduction channel by a distance that 
is large compared to the graphene lattice constant but 
small compared to the Fermi wave length in the regime 
studied experimentally (|n| < 7 x 10i2[cj^-2] |^^[^^ 
act like Coulomb scatterers in the MDF model but will 
not produce strong intervalley scattering. 

The difference between Coulomb and short-range cr™'" 
values can be rationalized by the following argument. In 
the short-range case, the golden-rule relaxation time for 
delta-function scatterers diverges (tq cx l/|i?_F| oo) 
at the Dirac point and kpi remains finite, which sug- 
gests that standard transport theory considerations may 
apply at least qualitatively once the non-vanishing den- 
sity of states in disordered systems is accounted for. In 
contrast, Boltzmann transport theory suggests that the 
Dirac point of the Coulomb scatterer model is in the 
strong disorder limit because Tc oc \Ep\ — s- 0. We ex- 
amined this idea by varying the carrier density n and the 
ionized impurity density Ui independently, finding that 
the conductivity for Coulomb scattering model appears 
to be a function of n/ui only. Letting n — > at finite 
Hi and letting Ui ~f oo at finite n are equivalent be- 
cause (7 (t""" ^ e'^/h. This finding is consistent with 
the idea that at the Dirac point MDF Coulomb scatter- 
ers are always in the strongly disordered limit. As long 
as these states remain delocalized, however, the conduc- 
tivity cannot vanish; the non-zero conductivity at the 
Dirac point is purely due to the quantum mechanical na- 



ture of delocalized states. Although this argument is not 
quantitative, it makes it clear that there is an essential 
difference between the Coulomb and short-range cases. 
A related difference is also seen in the inset in Fig. 1(a) 
which compares densities of states near the Dirac point 
{E = 0) for short-range scatterers (green dashed line) 
and Coulomb scatterers (blue solid line). The prominent 
dip at i? = in the short-range case is replaced by a 
smooth minimum at a larger value in the Coulomb case. 
One interpretation of the increase in Dirac point density 
of states is that the carrier density fluctuates spatially 
in the smooth Coulomb potential. Non-zero local car- 
rier densities at the Dirac point could explain both the 
increase in density of states and the increase in conduc- 
tivity. For short-range scatterers, our simulations have 
focused on the Boltzmann dominated kpS >> 1 regime. 
We note that systems with kpi ^ 1 oi smaller behave 
more like Coulomb scatterer systems at low densities, 
and have larger minimum conductivities and densities of 
states. 

Although our theory does not account for electron- 
electron interactions, the issue of screening effect near 
the Dirac point requires comment. We use the T = 
Thomas-Fermi approximation U{q) — 2Tre'^/e{q + 
Sttc^Pf) to give an indication of the importance of screen- 
ing, although this approximation must fail for Ep —^ 
since it applies strictly only for disorder potential range 
much longer than the Fermi wavelength and ksT <ti Ep. 
The density of states at the Dirac point pp vanishes in the 
clean limit whereas it can be evaluated self consistently 
in the disordered case, leading to the small but finite 
value seen in the inset in Fig. 1(a). The inset in Fig.l (c) 
compares the conductivity for unscreened and screened 
case at several different system sizes. We find that the 
conductivities evaluated with screened interactions are 
suppressed somewhat compared to the unscreened case. 
It appears likely that the precise value of cr™'" may de- 
pend indirectly on electron-electron interactions. 

We emphasize that we are able to evaluate the Kubo 
formula only over a relatively small range of finite sys- 
tem sizes. If we use the mobility to convert to physi- 
cal length units the minimum system size we study is 
~ 0.1/im, not much smaller than the size of the crystal- 
lites studied experimentally, [l^l Over the range of sizes 
we are able to study the conductivity does increase very 
slowly with system size, even at E = 0, as shown in 
the inset of Fig. 1(c). This weak size dependence might 
be hinting that the infinite conductivity predicted for 
two-dimensional systems with symplectic svmmetrv[29j| 
by scaling arguments will emerge in the MDF model at 
extremely large system sizes even at _E = 0. It is not 
obvious to us that the scaling theory applies at the Dirac 
point for either Coulomb or short-range scatterers since 
there is no length scale on which a ^ e'^/h. Alternately, 
the weak size dependence we find might be related to en- 
hancement of carrier-density spatial fluctuations, which 
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can be imagined in terms of puddles of electrons and 
holes, at larger system sizes in the long-range disorder 
potential case. Since Dirac fermions are not strongly 
localized, these puddles are not isolated but effectively 
percolate [3^ and are constrained by boundary condition 
in finite size systems. The large length scale limit of 
the conductivity at zero temperature is more obvious for 
real graphene than for its (two-component) MDF model, 
since intervalley scattering will always become relevant 
and lead to localization. [la. [itI ITst 

Finally we comment on statistical properties of the 
conductivity and the resistivity at the Dirac point. Fig- 
ure. 1 (d) shows the distribution function of the con- 
ductivity (main panel) and the resistivity (inset). The 
distribution function of a has a sharp peak near a ~ 
{l/'K)e'^ /h, and a large tail on the the large a regime. On 
the other hand, the resistivity distribution function has 
a broad peak near /i/e^, as seen in experiment . [2^ 
Discussion - We find that a MDF model with Coulomb 
scatterers is able to account for two key findings of exper- 
iments on graphene sheets, namely that the conductivity 
is proportional to carrier density away from the Dirac 
point and that the minimal conductivity per channel is 
finite and larger than the SCBA value obtained theo- 
retically for a MDF model with short-range scatterers. 
The impurity density rii in the Coulomb model should 
be associated with the density of ionized impurities that 
are located in the substrate within a Fermi wavelength 
of the graphene plane. In this interpretation, the mobil- 
ities measured in current samples [2^ [2^ correspond to 
rii ~ 5 X 10^^ [cm^]. The property that the conductiv- 
ity is proportional to carrier density in the Boltzmann 
regime suggests dominant smooth intravalley scattering 
that could be Coulomb in character. If so intervalley 
scattering is likely to be irrelevant in graphene in the 
Boltzmann transport regime, although it will always be 
important at finite carrier densities in the weak local- 
ization regime when the temperature is low enough that 
the phase coherence length exceeds the intervalley scat- 
terling length, Zq oc 1/|£'f|. When intervalley scattering 
is weak, the MDF model applies and momentum space 
Berry phases change constructive interference of back 
scattering particles into destructive interference 15. 3l| . 
implying weak antilocalization. When local scatterers 
dominate in graphene intervalley and intravally scatter- 
ing rates IJ, Il5l . Il7l | are comparable and the accumu- 
lated Berry phases are randomized , so that all states are 
localized. |l3|. l2Cll| Qualitative studies of electron-phonon 
interaction[32] and electron-electron interaction[lo| in 
graphene give important information on the coherence 
length. In a weak magnetic field, the magnetic length 
limits the coherence length. One recent magnetoresis- 
tance study has shown that highly doped graphene ex- 
hibits a negative magneto resistance indicating weak- 
localization, while states remain delocalized in the vicin- 
ity of the Dirac point [2^ . These properties further sup- 



port the view that graphene is described by a MDF model 
near the Dirac point. Our numerical results suggest that 
the MDF conductivity is finite and ^ jh at the Dirac 
point over the relevant range of the system size and the 
coherence length. 
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Note added ~ After submitting the present paper, 
we became aware of related work in which the impor- 
tance of the long-range nature of disorder scatterers is 
discussed. 331 
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